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Abstract 

Two new families of completely integrable perturbations of the A^-dimensional 
isotropic harmonic oscillator Hamiltonian are presented. Such perturbations de- 
pend on arbitrary functions and free parameters and their integrals of motion 
are explicitly constructed by making use of an underlying /ig-coalgebra symmetry. 
Several known integrable Hamiltonians in low dimensions are obtained as par- 
ticular specializations of the general results here presented. An alternative route 
for the integrability of all these systems is provided by a suitable canonical trans- 
formation which, in turn, opens the possibility of adding (N — 1) 'Rosochatius' 
terms that preserve the complete integrability of all these models. 
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1 Introduction 



The harmonic oscillator is indeed the paradigm of integrable Hamiltonian dynamics 
on the A^- dimensional (ND) Euclidean space, and its superintegrability properties are 
very well established through group-theoretic methods (see, for instance, [Ill2])- There- 
fore, the construction of non-trivial integrable perturbations/deformations of the ND 
harmonic oscillator (understood in a broad sense as the addition of any ND potential 
term -not necessarily small- that preserves the complete integrability of the system) 
are interesting from both physical and mathematical viewpoints. 

In this respect, although certain 2D and 3D perturbations have been thoroughly 
studied (see, for instance [3]-|ll] and references therein) the number of ND results is 
more limited. Among them we can quote the search for coupled polynomial pertur- 
bations [I^-llSj, Garnier systems [T^], the Smorodinski-Winternitz Hamiltonian ^201 
UT\ . Calogero-Moser systems [221 [23] and coupled oscillators coming from symmetric 
spaces [23]. From a more generic perspective, integrable oscillators on ND curved 
spaces can also be considered as deformations of the ND Euclidean ones in which the 
kinetic energy term is also perturbed through functions of the canonical coordinates 
containing a curvature parameter (see [25l|26l [27] and references therein). 

In this paper we demonstrate the complete integrability of the following ND natural 
Hamiltonian system 

j=l i=l \i=l i=l J \i=l J 

where V_2 is any homogeneous function of degree —2 in the canonical coordinates and 
such that V_2 can be written in terms of the variables J2iLi •^ili ^^"^ Yld=i if- Here, 
6i and Aj {i = 1, . . . , A^) are (A^ + 1) free parameters and J' is an arbitrary radial 
perturbation. Since the expressions for the integrals of the motion are analytic in 6i, 
the oscillator term can also be removed without altering the integrability properties 
of the system. As a consequence, by considering appropriate specializations of the 
dimension A^, the arbitrary functions V_2 and and the (A^ + 1) free parameters 
Xi {i = 1, . . . , N), we can obtain many different integrable systems. 

In Section 2 we shall prove this result by giving explicitly a set of (A^— 1) functionally 
independent integrals of the motion for "Hr, hat will be obtained by making use of the 
underlying /ig-Poisson coalgebra symmetry [281 (291 [30] of the Hamiltonian Tir- In 
fact, such symmetry implies that this Hamiltonian has two different sets of (A^ — 1) 
independent integrals [31], and among both sets we find a total number of (2A^ — 4) 
independent constants of the motion in involution with Tir, which implies that fll.ip is 
a superintegrable system (albeit not maximally superintegrable). 

Moreover, another family of completely integrable perturbations of the ND oscillator 
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is introduced in Section 3, namely 

'h = -J2p"+^^J2^"+j'[ E E^^^O (1-2) 

i=l i=l \l<i<j / \i=l J 

where J-" and Q are arbitrary functions and 6i and A, = 1, . . . , N) are again (A^+1) free 
parameters. Again, the second term {J-" + Q) can be interpreted as two independent 
famihes of integrable deformations/perturbations whose superposition preserves the 
complete integrability of the whole system. The (super) integrability of this family of 
ND Hamiltonian systems can also be proven by restorting to its /le-Poisson coalgebra 
symmetry and we shall present a number of known integrable polynomial perturbations 
that can be easily obtained as particular cases of H. 

We would like to mention that in the 3D case Hamiltonians f ll.ip and (11.21) are 
subclasses of the systems considered in [HI [ID] through a 3D realization of a so-called 
'extended' sZ(2,M) Poisson-Lie algebra (which is just isomorphic to h^). However, the 
ND integrable generalization here presented is only possible by endowing with the 
appropriate underlying coalgebra structure and by imposing the full Hamiltonian to 
be /ig-coalgebra invariant. 

In Section 4 we discuss how the complete integrability of (II. 2p can be alternatively 
explained through the introduction of a ND canonical transformation {q,p) — )■ {Q,P) 
under which the complete integrability of (II. 2p comes immediately from both its sepa- 
rability in terms of the (new) Qi coordinate and the radial symmetry of the remaining 
(A^ — l)-dimensional Hamiltonian. Moreover, we show that this radial symmetry can 
be broken without destroying the complete integrability of the system through the 

TV 

addition of a 'Rosochatius' contribution [321 ES] of the form ^ h/Qf- This fact will 

i=2 

be algebraically rephrased through the introduction of a new ND symplectic realiza- 
tion for the /ig-Poisson coalgebra that includes the previous Rosochatius terms. As a 
consequence, a further integrable generalization of the Hamiltonians (II. ip and (11.20 is 
presented. 

2 A new ND integrable Hamiltonian 

In order to prove the complete 
the functions 

N 

i=l 
N 

i=l 



integrability of the ND Hamiltonian (11.10 we consider 



AT ^ / ;^2\ 

A- = J2 ^ili ^ = I] ( - Y ) 

i=l i=l ^ ^ 

N N 

B- = Y.i M = Y.Xl (2.1) 



i=l i=l 
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where A, are arbitrary parameters and N fixes the number of degrees of freedom. These 
six functions provide a ND symplectic reahzation of the six- dimensional /le-Poisson 
coalgebra with Lie-Poisson brackets given by 



(2.2) 



{K,A+} = A+ {K,A^} = -A^ {A_,A+} = M 

{K,B+] = 2B+ {K,B^} = -2B_ {B_, B+] = AK + 2M 

{A+,i?_} = -2A_ {A+,fi+} = {M, ■} = 

{A_,B+] = 2A+ {A_,i?_} = 0. 

The Casimir functions for [35] are given by the central generator M and by 

C = {MB+ - Al){MB^ - Al) - {MK - A^A+ + M'^/2f. (2.3) 

Now, let 

(N ^ \ 

i=l i=l / 

be a homogeneous function of degree —2 in the canonical coordinates that can be 
written in terms of the generators A_ and Then the ND Hamiltonian 

N N / N ^ \ f ^ \ 

= ^ E + ^1 E ^ + ^-AY. E + ^ E ^0 ^^-^^ 

i=l i=l \i=l i=l / \i=l J 

where JT" is an arbitrary radial function and 5i, Aj {i = 1, . . . , A^) are free parameters, 
is completely integrable. 

In order to prove this result, we realize that Hr is /ig-coalgebra invariant, since it 
can be written as a smooth function of the /le-coalgebra generators in the form 

nr = \b+ + (5i B_ + V_2 (A-, B^) + J (B_) . (2.6) 

This invariance immediately implies (see [30]) that the following (A^ — 2) functionally 
independent functions 

m 

= E ^i^Pi^k - Pkqj)+>^j{pkqi - Piqk)+h{piqj - Pjqdf m = 3,4,. . . ,N 

i<j<k 

(2.7) 

Poisson-commute with T-Lr and among themselves. In fact, a given integral C[^^ is just 
the Casimir function (12. 3p divided by M and expressed in terms of the corresponding 
m-particle realization of the algebra given in the form ( 12. ip . 

Therefore, we have to find only one more independent integral in order to get 
complete integrability. In this case such an additional integral can be found through 
direct computation by making use of the Poisson brackets (12. 2p . and reads 

Ir = B_B+ ~{K+ ^Mf + 2B_ V_2 {A_, B_) (2.8) 
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that written in terms of the canonical coordinates gives: 

l<i<j \i=l ) \i=\ t=l / 

Again, the involutivity of "H^ and Xr with respect to the integrals (12 .yp is a direct 
consequence of the /le-coalgebra symmetry of both functions and the bracket {"Hr-jX^} = 
is easily proven by direct computation. Finally, the functional independence of can 
be easily proven in the harmonic oscillator case J = V_2 = with Aj = 6ij for a given 
j. Therefore, the functional independence follows for the generic case with arbitrary 
Xi and non-vanishing J' and V_2 functions. 

Moreover, the /le-coalgebra symmetry implies that Tir is not only integrable but 
superintegrable, since there exists an alternative set of integrals [30] 

N 

C/i6,(m) = {KiPjQk - PkQj) + XjiPkQi - PiQk) + XkiPiQj - PjQi)) (2.9) 

N-m+l<i<j<k 

where m = 3, . . . , N. Note that C/ig^(iv) = so Hr has a total number of {2N — 4) 
functions in involution (including X^). 

We recall that the ND isotropic oscillator (without any perturbation) has a total 
number of {2N — 2) independent integrals, that make it a maximally superintegrable 
system. From this perspective, the perturbation given by J-" and Q strongly preserves 
the integrability properties of the oscillator Hamiltonian, since only two integrals of 
the motion are lost under perturbation as a consequence of imposing the invariance 
of the system. 

Remark. Note that in the case = 3 the Hamiltonian (12. 5p is a particular subclass 
of the Hamiltonians presented in [9]. This is due to the fact that in our construction 
we have imposed V_2 to be /ig-coalgebra invariant, which is an additional constraint 
with respect to [9], where the only condition imposed on V_2 was to be a homogeneous 
function of degree —2 in the 3D canonical coordinates. 



2.1 A new ND integrable model 

If we take 



V,.(A-,fl-) = . ^_g^-/_-^,^ (2.10) 



and Ai = A2 = ■ • • = A V = 1 , we obtain the Hamiltonian 



N 

N N J2 QiQj / N 



«.-^Ep? + ^-E^? + ^ w + (2-11) 

\i=l J l<i<j 



^ i=l i=l I \^ n '^ \ \^ (n.-n.Y \ i=l 
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in which we see that the function V_2 (^-, gives rise to an homogeneous function of 
degree —2 in the canonical coordinates. In the N = 2 case (and with a vanishing radial 
perturbation J') the quantum mechanical analogue of (12.111) has been studied in [36] 
as an exactly solvable 'generalization' of the 2D Calogero-Moser model. Therefore, 
presumably (12. lip provides a good candidate for a new exactly solvable (classical and 
quantum) model in arbitrary dimension. 



3 Another ND integrable perturbation 

Many other possibilities for a perturbed ND oscillator Hamiltonian defined on the 
he Poisson coalgebra can be considered, but only very few special choices admit the 
existence of an additional integral X. This is the case of the function 

n = ^B+ + 6iB_+T{Cg_)+g{A_) (3.1) 

where the generators of are taken in the ND symplectic realization (12. ip and 

N 

Cg_ = MB^ -Al= i^jli - ^iljf (3-2) 

i<i<j 

is just the Casimir for a (1+1)-Galilei subalgebra of spanned by {B^, A^, M}. 
Therefore, we get the Hamiltonian 

^ = 2^^' + ^!^^' + -^ E (3.3) 

i=l i=l \l<i<j J \i=l J 

for which the additional integral X can also be found in an /ig-invariant form through 
direct computation and reads: 

I = Al + 26iA'i + 2g (A_) M. 

In terms of the ND symplectic realization, this last integral is explicitly written as 

(N \2 / N \2 / ^ \ / ^ \ 

A.pJ+2 5if ^ X,qA + 2 f ^(E^^^O • (3.4) 

Note that X does not depend on the choice of the function J-". 

We stress that the /le-coalgebra symmetry ensures that X will Poisson-commute 
with all the C[^^ and Cft6,(m) integrals, and the fact that {'H,X} = is easily proven by 
direct computation. The functional independence of X with respect to the remaining 
integrals is proven by taking 6i = and Q = 0. In this case X = A"^ only depends on 

the momenta, and is indeed functionally independent of C[^^ (that does not depend on 
5i, T and ^) and of the Hamiltonian. This functional independence will be mantained 
under the 'deformations' (5i 7^ and ^ 7^ 0. 
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3.1 Some examples 



The Hamiltonian ( 13.3P includes many different systems, and for = 2, 3 some of them 
can be found in the hterature. For instance, if we consider J-'(x) = a/x we get 

^ N N / ^ \ 

^ = ^Y.P^ + ^^Y.^" + — + G[J2X.qA (3.5) 

i<i<j 

which in the case = 2 is a rational Calogero-Moser system plus an arbitrary function 
Q. The integral of the motion in this case is f l3.4p with N = 2. Note that the integrals 
(12. 7p and (12. 9p only exist from = 3 and X does not depend on the choice of the 
potential J-". 

Nonlinear polynomial oscillators can be obtained from (I3.3P through a perturbation 

k I 

= J^{x) + g{y) = ^ a, + (3-6) 

i=2 j=3 

thus obtaining an integrable superposition of a polynomial deformation of degree 2k 
in the coordinates (the argument x = Cg_ is quadratic in q) and another one of degree 
I (the argument y = A- is linear in q). If we analyse separately the homogeneous 
perturbations with a fixed degree m, for odd m only g{A_) can contribute, which 
means that the richest cases will be obtained for m even. 

In particular, the homogeneous quartic case is obtained by taking cj) = a2 x"^ + (3^ y^. 
In the 2D case, if we take Ai = 1, A2 = we obtain a separable Hamiltonian. On the 
other hand, by taking Ai = 1, A2 = 1, 0:2 = /Ss = e/2 we recover the well-known 1:6:1 
integrable quartic perturbation (case (4)3' in [8] with = 0): 

n^'^ = \{pI+ pI) + 5i (g? + ql) +e{qt + qt + H^l) ■ (3-7) 

This system is known to be connected with the Hirota-Satsuma coupled KdV system 
(see |37]). In three dimensions, the case Ai = 1,A2 = 0, A3 = again provides the 
separability of the qi coordinate in the form 

7^(3) = i(p2 + p2 ^ ^ ^^^^2 ^ ^2 ^ ^2) ^ ^^^^2 ^ ^2)2 ^ ^^^4 (3 

which is also reflected in the invariant X = pl + 2((5ig^ + f^^qf). Moreover, the second 
invariant reads Cj^^ = {p2q3 — P3Q'2)^; which indicates the axial symmetry with respect 
to the coordinates {q2,q3)- This is the case (1) of Table 1 in [11]. On the other hand, 
if we take Ai = A2 = A3 = 1 we get 

+2 a2 ((92 - gi)' + (gs - qiY + (^3 - ^2)') + /^s (gi + g2 + ga)' (3.9) 
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which is just the 3D Chudnowski potential (see [9l[38]) plus an additional quartic term. 

The generic homogeneous sextic case is obtained by taking cf) = a^x^ + y^. In 2D, 
separability is trivially obtained when A2 = 0. The case Ai = A2 = 1 with a2 = = h/2 
gives the known coupled sextic oscillator (see [T7], case (i) of Table 5.3.): 

H^'^ = \{pI+ P^) + 5i [ql + ql) +h{ql + ql + 15 qlqM + ql)) . (3.10) 

In 3D the generic integrable homogeneous sextic perturbation coming from ( 13. 3p gives, 
under the choice Ai = 1, A2 = 0, A3 = 0, a system which is again separable in qi and 
whose axial symmetry in {qi,q2) is also recovered. As we shall see in the following 
section, this separability pattern can be explained in general for the Hamiltonian (13. 3p . 



4 A canonical transformation 



Without any loss of generality, let us assume that Ai 7^ and A^ 7^ for at least one 
k ^ 1. Let us consider the N x N orthogonal matrix R = (ciij) with matrix elements 

A, 
A 



\(^\ _\) ^^'"^ '^■> ■ ■ ■ ■> 



N 

where A = 1 / ^ A^ and Sij is the Kronecker delta. Then R defines a canonical transfor- 

i=l 

mation on the (column) vectors q := (gi, q2, . . ■ , ^at) and p := {pi,P2, ■ ■ ■ ,Pn) through 

Q = R.q, P = R-p (4.11) 

where 

i?* = R-^ = R, detR = -1 (4.12) 

and the (column) vectors Q := (Qi, Q2, ■ ■ ■ ,Qn) and P := (Pi, P2, . . . , Pn) are the new 
canonical variables fulfilling {Qi,Pj} = 6ij. 

In terms of these new variables the A^-particle symplectic realization of reads 

A^ 



A+ = XPi A_ = XQi K = Y^QiPi 

N N 

B^ = Y.P' B_ = Y,Ql M = \\ (4.13) 

i=l i=l 

and the integrals of the motion coming from the Casimir are transformed into 

m 

CM^X' J2 {PjQ, - PiQjf m = 3,...,N. (4.14) 

2<i<j 



Essentially, the existence of this canonical transformation shows that the symplectic 
realization of (12. ip with A2 = A3 = ■ ■ ■ = Aat = is the generic one. In particular, 
under this canonical transformation the Hamiltonian (13. 3 p reads 

^ N N / ^ \ 

^' = ^ E + E + E + ^ (4.15) 

i=l i=l \ i=2 / 

whilst the additional integral (13. 4p is just 

I' = X^P^ + 26iQl+2g{XQi)}. (4.16) 

Therefore, the {N — 2) integrals (14.141) just reflect the radial symmetry of the Hamilto- 
nian T-L' with respect to the new variables i = 2, . . . , N, and the remaining one (14.161) 
expresses the separability of the i = 1 degree of freedom. 

At this point it worth stressing that the canonical transformation R allows us to 
deduce immediately -and without making use of any invariance- the existence of 
the invariants Lij = {PjQi — PiQj) with 2 < i < j < N, which are linear in the 
momenta. Nevertheless, such linear integrals are not generically in involution, and 
the construction provides automatically the subfamily (14.141) of (A^ — 2) quadratic 
integrals in involution. Anyhow, the linear 'angular momentum components' L^j are 
indeed useful in order to simplify the study of the dynamics of the system and to 
construct its explicit solutions. 

Under the canonical transformation R, the Hamiltonian (11.11) is also mapped into 

^: = ^ E + ^1 E + ^ E + ^ E (4-17) 

i=l i=l V «=2 / \ t=l / 

and its additional integral is written as 

x:.= 5^(p,g.-p.Q,)2+2 K^Qnv_2 Agl,^Q^ (4.18) 

l<i<j \i=l ) \ i=l / 

Note that in general (I4.17P is no longer a separable system, although it does exhibit 
radial symmetry with respect to the i = 2, 3. . . . , A^ degrees of freedom. 

4.1 Rosochatius terms through a new symplectic reaUzation 

After the previous discussion it comes out that the complete integrability of the Hamil- 
tonian (13.31) is evident when written in the new variables {Q,P). Nevertheless, this 
result opens the path to the construction of a non-trivial completely integrable general- 
ization of both (I4.15P and (I4.17p . since if radial symmetry appears then a Rosochatius- 
type potential can be added without altering the integrability properties of the Hamil- 
tonian (see [25j, [27j and [3lj for a detailed discussion in terms of the appropriate 
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symplectic realization of the sl{2) Poisson-coalgebra) . This enables us to consider a 
modification of the iV-particle symplectic realization ( I4.13P including a set of (A^ — 1) 
dimensional Rosochatius terms within the S+ generator as follows: 

A+ = APi A^ = XQi K = Y,QiPi-^ 

i=l 

Note that when bi > the Rosochatius terms can be properly considered as centrifugal 
ones. It is immediate to prove that fl4.19p is a symplectic realization for h^, and the 
Casimir functions under this realization read 



[ 1=2 2<i<j 2<i<j ^ Vi / J 



m = 3,...,N. 

(4.20) 



If we now use this symplectic realization for the Hamiltonian fl3.ip we get 

^ N N / ^ \ ^ /i 

^' = ^E^'+^iE^'+'^h'E^n+^(^<5i) + E^ (4.21) 

i=l 1=1 \ i=2 J i=2 ^* 

which is again completely integrable with constants of the motion given by fl4.16p and 
( I4.20p . Note that if the canonical transformation is reversed, the latter Rosochatius 
terms would be written in a non-trivial way in terms of the original coordinates 
(gi,g2, • • • ,gAr)- Namely, 



Qi = ^\y M^\ (4-22) 



Q\ = y^2^y^^_ y^^y [ ' Ai)(A, gi + A q^) - A. ^ \,q, ) 2 > 2. (4.23) 



In the same way, the symplectic realization fl4.19p can be substituted onto the 
'abstract' Hamiltonian H!^ (12.61) and leads to 

1^ Af / ^ \ ( ^ \ ^ h 

= ^ E + ^1 E + ^-M S + ^ + 5Z ^- (4-24) 

i=\ i=\ \ i=\ J \i=l J 1=2 

This is indeed a new completely integrable Hamiltonian with constants of the motion 
in involution given by (I4.20p and by the symplectic realization (I4.19p of 

Xr = B_B+ -{K+ ^Mf + 2B_ V_2 {A_, B_) (4.25) 
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that written in terms of the new canonical coordinates reads: 



l<i<j 1=1 j=2 \i=l / \ 1=1 

AN) jv 



2 



A2 



thus showing its functional independence with respect to the integrals C^™"* f l4.20p . 

Summarizing, we have presented a constructive approach of two families of com- 
pletely integrable perturbations-deformations of the ND harmonic oscillator. The fact 
that all these Hamiltonians can be written as functions of the /ig-Poisson coalgebra 
under suitable iV-particle symplectic realizations makes it possible to demonstrate the 
(super) integrability properties of both classes of systems. Among all the Hamiltonians 
here presented, the model ( 12. lip seems to deserve some attention. In particular, its 
explicit solvability could be addressed through the 'cluster variables' technique [2S] in 
order to select the appropriate dynamical variables by making use of the underlying 
coalgebra structure. Another interesting possibility for further work would be the con- 
struction of the analogues on ND spaces with constant curvature of the Hamiltonians 
here presented. This could be feasible by mimicking the s/(2,M) approach described 
in P5| 1221 but now making use of the corresponding /ig-coalgebra formalism. 
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